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Abstract 

In this paper we consider the relation between symmetries and first integrals for 
both continuous canonical Hamiltonian equations and discrete Hamiltonian equa- 
tions. We observe that canonical Hamiltonian equations can be obtained by vari- 
ational principle from an action functional and consider invariance properties of 
this functional as it is done in Lagrangian formalism. We rewrite the well-known 
Noether's identity in terms of the Hamiltonian function and symmetry operators. 
This approach, based on symmetries of the Hamiltonian action, provides a simple 
and clear way to construct first integrals of Hamiltonian equations without integra- 
tion. A discrete analog of this identity is developed. It leads to a relation between 
symmetries and first integrals for discrete Hamiltonian equations that can be used 
to conserve structural properties of Hamiltonian equations in numerical implemen- 
tation. The results are illustrated by a number of examples for both continuous and 
discrete Hamiltonian equations. 



1 Introduction 

It has been known since E. Noether's fundamental work that conservation laws of differential 
equations are connected with their symmetry properties [1]. For convenience we present here 
some well-known results (see also, for example, [21E1II]) for both Lagrangian and Hamiltonian 
approaches to conservation laws (first integrals). 
Let us consider the functional 

h(u) = / L(x,u,u\)dx, (1.1) 
Jn 

where m ) are independent variables, u = (it 1 , u 2 , u n ) are dependent variables, 

Ui = (uf) are all first order derivatives u\ = f2 is a domain in R m and L(x, w, u\) is a first 
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order Lagrangian. The functional (11.11) achieves its extremal values when u(x) satisfies the 
Euler-Lagrange equations 



6L dL ( dL , 

A hrr =°. k = i,...,n, (1.2) 



5u k du k \ du 
where 

D ' = w +vt w + <^ + -- i = 1 "" m 

are total differentiation operators with respect to independent variables x % . Here and below we 
assume summation over repeated indexes. Note that equations (11.21) are second order PDEs. 
We consider a Lie point transformation group G generated by the infinitesimal operator 

X = ?{ X ,u)^ + Tf(x,u)±; + ..., (1.3) 

where dots mean an appropriate prolongation of the operator on partial derivatives [5j El [7J [8] . 
The group G is called a variational symmetry of the functional L(iz) if and only if the Lagrangian 
satisfies pQ 

X(L) + LA(f) = 0, (1.4) 

where X is the first prolongation, i.e. the prolongation of the vector field X on the first 
derivatives u k . We will actually need a weaker invariance condition than given by Eq. (jl.4p . 
The vector field X is a divergence symmetry of the functional h(u) if there exist functions 
V l (x, u, Mi), i — 1, m such that [9] (see also [HI [3, E]) 

x(L) + LA(e) = A(n- (i.5) 

An important result for us is the following: If X is a variational symmetry of the functional 
L(m), it is also a symmetry of the corresponding Euler-Lagrange equation. The symmetry group 
of Eqs. (II. 2p can of course be larger than the group generated by variational and divergence 
symmetries of the Lagrangian. 

Noether's theorem [1] states that for a Lagrangian satisfying the condition (11 .4j) there exists 
a conservation law of the Euler-Lagrange equations (II. 2p : 

Di(?L+tf-evi)^=0. (1.6) 

This result can be generalized: If X is a divergence symmetry of the functional L(u), i.e. 
equation (II. 5p is satisfied, then there exists a conservation law 

a(^+(^-^«J)^-V*)=0 (1.7) 

of the corresponding Euler-Lagrange equations. 

The strong version of the Noether's theorem [7] states that there exists a conservation law 
of the Euler-Lagrange Eqs. (II. 2p in the form (II .6p if and only if the condition (11.41) is satisfied 
on the solutions of Eqs. (II. 2p . 
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In the present paper we are interested in canonical Hamiltonian equations 

q= — , Pl = ~ W . = l,...,n. (1.8) 
These equations can be obtained by the variational principle from the action functional 

5 [' ( Pi q i -H{t,q,p))dt = (1.9) 
Jti 

in the phase space (q, p), where q = (q 1 , q 2 , q n ), p = (p\,p 2 , ...,p n ) (see, for example, [TUl 

EU). 

Let us note that the canonical Hamiltonian equations (11. 8ft can be obtained by action of 
the variational operators 

S d d 

bp,, dpi dpi 

A- = JL- D JL ] i = l,..., n , (1.11) 

oq l dq l dq L 

where D is the operator of total differentiation with respect to time 



d ■ d d 

D =m + iW +pi W i + -' (L12) 

on the function 

p { q l - H(t,q,p). 

The Legendre transformation relates Hamiltonian and Lagrange functions 

L(t,q,q)= Pl <f-#(t,q,p), (1.13) 

where p = q = It makes it possible to establish the equivalence of the Euler-Lagrange 
and Hamiltonian equations @]. Indeed, from Euler-Lagrange equations for one independent 
variable (m = 1) 

5L dL fdL\ 

we can obtain the canonical Hamiltonian equations (11.81) using the Legendre transformation. 
It should be noticed that the Legendre transformation is not a point transformation. Hence, 
there is no conservation of Lie group properties of the corresponding Euler-Lagrange equations 
and Hamiltonian equations within the class of point transformations. 

Lie point symmetries in the space (t, q, p) are generated by operators of the form 

d d d 

X = at, ^P)g- t + V% q, P)^- + Ci(t, q, P)^-- (1-15) 

Standard approach to symmetry properties of the Hamiltonian equations is to consider so 
called Hamiltonian symmetries [6]. In the case of canonical Hamiltonian equations these are 
the evolutionary (£ = 0) symmetries (I1.15P 

X = n% q,p)^ + Ci(t, q, p)^ (1.16) 
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with 

l= dl_ c = i = l 

^ dpi ' dq 1 ' ' ' 

for some function I(t, q, p), namely, symmetries of the form 

1 dpi dq l dq % dpi 

These symmetries are restricted to the phase space (q, p) and are generated by the function 
I = I(t, q, p). For symmetry (11.181) the independent variable t is invariant and plays a role of 
a parameter. 

Noether's theorem (Theorem 6.33 in [6]) relates Hamiltonian symmetries of the Hamilto- 
nian equations with their first integrals. Restricted to the case of the canonical Hamiltonian 
equations it can be formulated as follows: 

Proposition 1.1 An evolutionary vector field X of the form U.16}) generates a Hamiltonian 
symmetry group of the canonical Hamiltonian system U.8\) if and only if there exists a first 
integral J(t,q, p) so that X = Xj is the corresponding Hamiltonian vector field. Another 
function I(t, q, p) determines the same Hamiltonian symmetry if and only if I = I + F(t) for 
some time-dependent function F(t). 

Thus, we obtain that the Hamiltonian symmetry determines a first integral of the canonical 
Hamiltonian equations up to some time-dependent function, which can be found with the help 
of these equations. This approach has two disadvantages. First, some transformations loose 
their geometrical sense if considered in evolutionary form (11.181) . Second, there is a necessity 
of integration to find first integrals with the help of ( 11. 17ft . In this approach it is also not clear 
why some point symmetries of Hamiltonian equations yield integrals, while others do not. 

In the present paper we will consider symmetries of the general form (11.151) . which are 
not restricted to the phase space and can also transform t. In contrast to the Hamiltonian 
symmetries in the form (11.181) the underlying symmetries have a clear geometric sense in finite 
space and do not require integration to find first integrals. We will provide a Hamiltonian 
version of the Noether's theorem (in the strong formulation) based on a newly established 
Hamiltonian identity, which is an analog of well-known Noether's identity for the Lagrangian 
approach. The Hamiltonian identity links directly an invariant Hamiltonian function with first 
integrals of the canonical Hamiltonian equations. This approach provides a simple and clear 
way to construct first integrals by means of merely algebraic manipulations with symmetries 
of the action functional. The approach will be illustrated on a number of examples, including 
equations of the three-dimensional Kepler motion. 

The paper is organized as follows: In section [2] we introduce a definition of an invariant 
Hamiltonian and establish the necessary and sufficient condition for H to be invariant. Sec- 
tion [3] contains the main propositions of present paper: Lemma 13.11 introduces a new identity, 
which is used in Theorem 13.21 to formulate the necessary and sufficient condition for existence 
of first integrals of Hamiltonian equations (Hamiltonian version of Noether's theorem in the 
strong formulation). In section HJ Lemma [4.11 introduces two more identities, which are used 
in Theorem 14.41 to formulate necessary and sufficient conditions for the canonical Hamiltonian 
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(1.17) 



(1.18) 



equations to be invariant. Section [5] contains example ODEs which are considered as both 
Euler-Lagrange equations and canonical Hamiltonian equations. In particular, we consider the 
equations of Kepler motion. In section [6] we present discrete Hamiltonian equations. Their sym- 
metries and first integrals are shown to be related in the same way as those for the continuous 
canonical Hamiltonian equations. Final section [7J contains concluding remarks. 



2 Invariance of elementary Hamiltonian action 

As an analog of the Lagrangian elementary action [El [7] we consider the Hamiltonian elementary 
action 

Pidq* - Hdt, (2.1) 
which can be invariant or not with respect to a group generated by an operator of the form fl 1 . 1 5 [) . 

Definition 2.1 We call a Hamiltonian function invariant with respect to a symmetry opera- 
tor U.15\) if the elementary action A2.1\) is an invariant of the group generated by this operator. 

Theorem 2.2 A Hamiltonian is invariant with respect to a group generated by the opera- 
tor U.15\) if and only if the following condition holds 

Qq l + PiD{rf) - X(H) - HD(£) = 0. (2.2) 

Proof. The invariance condition follows directly from the action of the operator X prolonged 
on the differentials dt and dq l , i = 1, ...,n [7j: 

X = {(*, q, p)| + n% q, P) A + Ut, q. p) A + fl®*^ + W*^. (2.3) 

Application of (12.31) to the Hamiltonian elementary action (12.11) yields 

X { Pl dq l - Hdt) = (Ci<f + PiD{rf) - X(H) - HD(0) dt = 0. 

□ 

Remark 2.3 From the relation 

L(t, q, q)dt = Pidq 1 - H(t, q, p)dt (2.4) 

it follows that if a Lagrangian is invariant with respect to a group of Lie point transformations, 
then the Hamiltonian is also invariant with respect to the same group (of point transformations) . 
The converse statement is false. For example, symmetries providing components of Runge- 
Lenz vector as first integrals of Kepler motion are point symmetries in Hamiltonian framework 
(point Hhty) . However, they are generalized symmetries in Lagrangian framework JE/. 

The proof follows from the action of operator (11.151) on relation (12.41) . 

Remark 2.4 The operator of total differentiation U.1S\) applied to Hamiltonian H and con- 
sidered on the solutions of Hamiltonian equations U.8\) coincides with partial differentiation 
with respect to time: 



dH .dH OH 

iTr + q ttj + Pi t; 



dH , , 

(2.5) 
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3 The Hamiltonian identity and Noether— type theorem 

Now we can relate conservation properties of the canonical Hamiltonian equations to the in- 
variance of the Hamiltonian function. 

Lemma 3.1 The identity 

(3.1) 



C4 l + PiD{rf) - X(H) - HD(0 = £ [ D(H) - ^ 



is true for any smooth function H = H(t, q, p). 

Proof. The identity can be established by direct calculation. □ 

We call this identity the Hamiltonian identity. This identity makes it possible to develop 
the following result. 

Theorem 3.2 The canonical Hamiltonian equations U.S\) possess a first integral of the form 

l = Vl rf-iH (3.2) 

if and only if the Hamiltonian function is invariant with respect to the operator M.15\) on the 
solutions of the equations U.8\) . 

Proof. The result follows from the identity (13. ip . □ 

Theorem 13.21 corresponds to the strong version of the Noether theorem (i.e. necessary and 
sufficient condition) for invariant Lagrangians and Euler-Lagrange equations [7]. 

Remark 3.3 Theorem \3.2\ can be generalized on the case of the divergence invariance of the 
Hamiltonian action 

+ p l D(r l l ) - X{H) - HD(0 = D(V), (3.3) 

where V = V(t, q, p). If this condition holds on the solutions of the canonical Hamiltonian 
equations U.8\) . then there is a first integral 

I = Pi r? -ZH-V. (3.4) 



4 Invariance of canonical Hamiltonian equations 

In the Lagrangian framework, the variational principle provides us with the Euler-Lagrange 
equations. It is known that the invariance of the Euler-Lagrange equations follows from the 
invariance of the action integral. The following Lemma 14.11 and Theorem 14.21 establish the 
sufficient conditions for canonical Hamiltonian equations to be invariant. 



6 



Lemma 4.1 The following identities are true for any smooth function H = H(t, q, p): 



^- (Ci<f + p,D(rf) - X(H) - HD(Q) 
DW)-?D { Q-x( a JL\ + ^(D ( H)- d 4\ (4.1! 



dpj J dpj \ dt J 



dpj \ dq l ) \ dpj J \ dp 



±(Q q * + Pl DW)-X(H)-HD(0) 



where the notation 5^ stands for the Kronecker symbol. 



Proof. The identities can be easily obtained by direct computation. □ 

Theorem 4.2 //a Hamiltonian is invariant with respect to the symmetry U.15\) . then the 
canonical Hamiltonian equations U.8\) are also invariant. 

Proof. For invariance of the canonical Hamiltonian equations fll.Sp we need the equations 

Dtf) - q j D(0 = X {^j , j = l,...,n 

D{d) - PjD^) = -X J = l,-..,n 

to hold on the solutions of the Hamiltonian equations [6]. These conditions follow from the 
identities (ED) and (Oji . □ 

Remark 4.3 The statement of Theorem \4 . S\ remains valid if we consider divergence symmetries 
of the Hamiltonian, i.e. condition because the term D(V) belongs to the kernel of the 

variational operators U.10\) . Fl.ll\) . 

The invariance of the Hamiltonian is a sufficient condition for the canonical Hamiltonian 
equations to be invariant. The symmetry group of the canonical Hamiltonian equations can 
of course be larger than that of the Hamiltonian. The following Theorem 14.41 establishes the 
necessary and sufficient conditions for canonical Hamiltonian equations to be invariant. 
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Theorem 4.4 Canonical Hamiltonian equations U.8\) are invariant with respect to the sym- 
metry U.15\) if and only if the following conditions are true ( on the solutions of the canonical 
Hamiltonian equations): 



5 

S Pj 



(( iq l + Pl DW)-X(H)-HD(0) 



JL(Qq l + ViDW)-X{H)-HD{i)) 



q=# p , p=-i? q 



q=H p , p=-Hq 



0, J = l,...,n, (4.3) 

0, j = l,...,n. (4.4) 

Proof. The statement follows from the identities (14. ip and (14.21) . □ 

It should be noted that conditions fl4.3[) and (14.41) are true for all symmetries of canonical 
Hamiltonian equations. But not all of those symmetries yield the "variational integral" of these 
conditions, i.e. 

(fcj* + Pi D{rf) - X{H) - HD{(-))\ . =Hp> . = _ H ^ = 0, 

which gives first integrals in accordance with Theorem 13.21 That is why not all symmetries of 
the canonical Hamiltonian equations provide first integrals. In the next section we illustrate the 
theorems, given above, on a number of examples. 



5 Applications 

In this section we provide examples how to find first integrals with the help of symmetries. 



5.1 A scalar ODE. 

As the first example we consider the second-order ODE 



« = 4 (5-1) 



U 



which admits Lie algebra L3 with basis operators 



Xi = L X 2 = 2t^- + u^, X 3 = t 2 ^- + tu-^. (5.2) 
at at ou at ou 



5.1.1 Lagrangian approach 

The Lagrangian function 



1 ( .0 1 



L{t J u,u) = -\u z --) } (5.3) 



which provides equation (15. ip as its Euler-Lagrange equation, is invariant with respect to X\ 
and X2. Therefore, by means of Noether's theorem there exist first integrals 

Jl = -i( i2 + i)' A = ««-<(" 2 + ^Y (5-4) 



S 



The action of the third operator X3 yields the divergence invariance condition 

XL + LD{£) = uii = D V (5.5) 
Due to the divergence invariance of the Lagrangian we can find the following first integral 

J 3 = -^ + (u-tu) 2 ). (5.6) 
Alternatively, one can find the last integral from another Lagrangian function 

'It \ 2 1 

t J U 



L(t,u,ii) = ( 7 -it) -3, (5.7) 



which is exactly invariant with respect to X%. 

It should be mentioned that independence of first integrals obtained with the help of the 
Noether theorem is guarantied only in the case when there is one Lagrangian which is invariant 
with respect to all symmetries. This condition is broken in the considered example. Therefore, 
the integrals obtained are not independent. Integrals (15.41) . (I5.6P are connected by the relations 



4J 1 J 3 -J 2 ' = 1. (5.8) 

Thus, any two integrals among (15.41) . (15.61) are independent. Putting J\ = A/ 2, J 2 = B and 
excluding ii, we find the general solution of the equations (15.11) as 

Au 2 + (At -B) 2 + 1 = 0. (5.9) 
5.1.2 Hamiltonian framework 

Let us transfer the preceding example into the Hamiltonian framework. We change variables 



8L 

q = u, p = — = u. 

ou 

The corresponding Hamiltonian is 



dL 1 / 9 1 \ . 

H &™) = *to- L= 2\f Jr ?)- (5 ' 10) 

The Hamiltonian equations 

1 

admit symmetries 

9 „ d d d o d d , .8 

Xi=«7, X 2 = 2t- + g— -p— , X 3 = t 2 - + tq— + (q-tp — . 5.12 
8t at 8q 8p at 8q 8p 

We check invariance of H in accordance with Theorem 12.21 and find that condition (12.21) is 
satisfied for the operators X\ and X 2 . Using Theorem 13.21 we calculate the corresponding first 
integrals 

1/0 i\ . . / o r 



h = -H = I p 2 + -A , h=pq-t[p 2 + - 2 ). (5.13) 
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For the third symmetry operator the Hamiltonian is divergence invariant with V3 = g 2 /2. In 



accordance with Remark I3.3[ it yields the following conserved quantity 

1 (t 

2 U 



h = — [- 2 + {q-tp) 2 ). (5.14) 



Note that no integration is needed to provide solutions of (15.111) . As we noted before in the 
Lagrangian case only two first integrals are functionally independent. Putting I\ = A/2 and 
I2 = B, we find the solution (15.111) as 

B — At 

Aq 2 + {At- Bf + 1 = 0, p = . (5.15) 

Q 

5.2 Repulsive one-dimensional motion. 

As another example of an ODE we consider one-dimensional motion in the Coulomb field (the 
case of a repulsive force): 

u=\, (5.16) 
which admits Lie algebra L 2 with basis operators 

* = 5- A ' 2 = 3 4 + 2 "!r (5 - 17) 



5.2.1 Lagrangian approach 

The Lagrangian function 



u 2 1 

L(t,u,u) = --- (5.1* 
2 u 



is invariant only with respect to X\. Therefore, Noether's theorem yields the only first integral 

u 2 1 

Ji = \ + ~- (5.19) 
2 u 

In this case the Euler-Lagrange equation admits two symmetries while the Lagrangian is in- 
variant with respect to one symmetry operator only. 

5.2.2 Hamiltonian framework 

We change variables 

dL . 



q = u, p= -wr = u 
ou 



and find the Hamiltonian function 



dL v 2 1 

H(t,q,p) = u—-L = + (5.20) 



The Hamiltonian equations have the form 



q = p, P = ~2- (5-21) 
10 



We rewrite symmetries in the canonical variables as the following algebra L2: 

d d d d 

Xt = - X 2 = 3t- + 2q— - p—. 
at at oq op 



(5.22) 



The invariance of Hamiltonian condition (12.21) is satisfied for operator Xi only. Applying 

(5.23) 



Theorem 13. 2\ we calculate the corresponding first integral 

/, = -* = -(£ + i 

V 2 q 

Application of operator X2 to the Hamiltonian action gives 

(q + P D( V )-X(H)-HD(Z)= P q 

Meanwhile, in accordance with Theorem 14.41 we have 
5 



2 q 



(5.24) 



dp 



((q + pD( V )-X(H)-HD(0) 



(q-p)\ 



-P+ -o 



(Cq+pD(rj)-X(H)-HD(0) 

1=P< V=^i 

We will show below that there exists a second integral of non-local character. 



It was shown in [12] that Eq. (15.161) can be linearized by a contact transformation. For 
equations (15.211) this transformation is the following 



p{t) = P{s) 



Q 2 (s) 



?(*)' 



4 , 



The new Hamiltonian 



H{s,Q,P) = -(P 2 + Q 2 



corresponds to the linear equations 



dQ =p ^ 
ds 



as 



(5.25) 
(5.26) 

(5.27) 



which describe the one-dimensional harmonic oscillator. These equations have two independent 
first integrals 



h = - 2 (P 2 + Q 2 ) 



+ s, 



I2 = arctan 

which let us write down the general solution of the equations (15.271) as 

Q = A sin s + B cos s, P = A cos s — B sin s, 



(5.28) 



(5.29) 



where A and B are arbitrary constants. Applying the transformation (15.251) to integral I2 we 
find the non-local integral for Eqs. (15.211) 



arctan 



( \ _ — 

\V2J V2~J to q 



dt 



7 3/2- 



(5.30) 
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5.3 Kepler motion 

Keplers problem is a special case of the two-body problem, in which the two bodies interact by a 
central force that varies in strength as the inverse square of the distance between them [21 [3J, S] - 
The three-dimensional Kepler motion of a body in Newton's gravitational field is given by the 
equations 

K 2 

q = p, p= — 3- q, r=|q|, q,peR 3 , (5.31) 

where K is a constant, with the initial data 

q(0) = q , p(0) = p . 
These equations are Hamiltonian. They are defined by the Hamiltonian function 

H( q ,p) = hp\ 2 -— . (5.32) 
2 r 

Among symmetries admitted by the equations (I5.3ip there are 

d Odd 
X = - X x = 3t- + 2q l — - Pl — , 
at at aq l dpi 

X — q j ® + q l ® p- ® +p- ® i^j 
lj ^ dq l ^ dqi ^ j dpi ^ 9pj' ^ 

d 

Yi = (2q l p k - q k pi - (q, p)*lfc)^ 

/ x 2 \ d 

+ ( piPk - (p, p)Sik - -^-(q q - (q, q)fo) J I = h 2, 3, 

where (f , g) = f T g is scalar product of vectors. 

The Hamiltonian function (15.321) is invariant for symmetries Xq and X^. Theorem 13.21 
makes it possible to find the first integral for symmetry X 

h = -H, 

which represents the conservation of the energy in Kepler motion. For symmetries Xij we 
obtain the first integrals 

hj = q l Pj ~ q J Ph i ^ j, 
which are components of the angular momentum 

L(q,p) = qxp. (5.33) 

Conservation of the angular momentum shows that the orbit of motion of a body lies in a 
fixed plane perpendicular to the constant vector L. It also follows that in this plane the 
position vector q sweeps out equal areas in equal time intervals, so that the sectorial velocity 
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is constant [I]. Therefore, Kepler's second law can be considered as a geometric reformulation 
of the conservation of angular momentum. 

The scaling symmetry X\ is not a Noether symmetry (neither variational, nor divergence 
symmetry) and does not lead to a conserved quantity. 

For each of symmetries Y\ the Hamiltonian is divergence invariant with functions 



Vi = q l ^(p, p) + ^-j - Pi(q, p) 



Hence, the operators Y\ yield the first integrals 



h = Q 1 ((p, p) - - Pj(q, p) 



/ = 1,2,3. 



1 = 1,2,3, 



which are components of the Runge-Lenz vector 

K 2 ( 1 \ 

A(q, p) = p x L - — q = q f #(q, p) + -|p| 2 J - p(q, p). (5.34) 

Physically, vector A points along the major axis of the conic section determined by the orbit 
of the body. Its magnitude determines the eccentricity [13]. 

Let us note that not all first integrals are independent. There are two relations between 
them given by the equations 

A 2 - 2#L 2 = and (A, L) = 0. 

The two-dimensional Kepler motion can be considered in a similar way. Let us remark 
that symmetries and first integrals of the two-dimensional Kepler motion can be obtained by 
restricting the symmetries and first integrals of the three-dimensional Kepler motion to the 
space (t,q 1 ,q 2 ,px,P2)- As the conserved quantities of the two-dimensional Kepler motion one 
obtains the energy 

1 K 2 

H(q,p) = -\p\ 2 - — , r=|q|, q=(<?\g 2 ), P = (Pi,p 2 ), 
one component of the angular momentum 

L 3 = q 1 P2 - q 2 Pi 

and two components of the Runge-Lenz vector 

M = q 1 ( H {q., p) + ^IpI 2 ^) - Pi(q, p)- 

A 2 = q 2 ^ff(q, p) + i|p| 2 ^ - p 2 (q, p). 

There is one relation between these conserved quantities, namely 

A\ + A\ - mh\ = K A . 

Further restriction to the one-dimensional Kepler motion leaves only one first integral, which 
is the Hamiltonian function. 
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6 First integrals of discrete Hamiltonian equations 



It is known that the preservation of first integrals (conservation laws) in numerical work is of 
great importance (see, for example, [MJC15]). Therefore, it makes sense to establish a discrete 
analog of the results presented for the continuous Hamiltonian equations. An analogous dis- 
crete framework would allow one to construct numerical schemes with first integrals 
for various applied problems. 



6.1 The discrete version of Hamiltonian action 

We will consider finite-difference equations and discrete Hamiltonians at some point (t, q, p) 
of a lattice. Generally, the lattice in not regular. The notations are clear from the following 
picture: 




To consider discrete equations we will need three points of a lattice. Prolongation of Lie 
group operator (11. 15ft for neighboring points (t_,q_,p~) and (t+,q + ,p + ) is the following: 



dt dq i l dpi ~ dt- ~ dql_ 1 dp~ 



(6.i; 



where 



£- = f(*-,q-:P~ 
= £0+,q+,p H 



77i = 77 i (t_,q_,p ), Q =C*(*-.q-,p" 
r)\ = 7]\t+, q+,p + ), c + = C(^+,q+,p H 



Hamiltonian equations can be obtained by the variational principle from the finite-difference 
functional 

(6.2) 



= $>^+ - - n(t,t + ,q, P + )h + ). 



Indeed, a variation of this functional along a curve q % = <pi{t), Pi = ipi{t), i = 1, ...,n at some 
point (£, q, p) will effect only two terms of the sum (16.21) : 

U h = ...+pM l - ql) ~ W(t-,t,q-,p)A_ +p+(q\ - q z ) - H(t, t+, q, p + )h + + ... (6.3) 
Therefore, we get the following expression for the variation 

SEh = 6 -^S Pi + S -P-6q l + *£st, (6.4) 
dpi bq % dt 
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where 5q l = fySt, 5pi = fySt, i = 1, n and 



(6.5) 



where ft = ft(t, £ + , q, p + ) and ft~ = H(t_,t, q , p). 

For the stationary value of the finite-difference functional (16. 2p we obtain the system of 
2n + 1 equations 

5H bU . 5H 

^ = °' *? = ' « = 1 '->», ^ = °- 

Thus, we arrive at the system of 2n + 1 equations 



D(q t ) = 1 - T , D(Pi) = -- ST , 1 = 1,... ,n, 



(6.6) 



which we will call discrete Hamiltonian equations. For convenience we use the following total 
shift (left and right) operators and corresponding discrete differentiation operators: 

S - 1 

Sf(t) = f(t ± ), D = ^—. 
±h ±h ±n± 

Let us note that the first 2n equations (16.61) are first-order discrete equations, which correspond 
to the canonical Hamiltonian equations (ll.Sp in the continuous limit. The last equation is of 
a second-order. Its continuous counterpart (see Remark 2.4) is automatically satisfied on the 
solutions of canonical Hamiltonian equations. In discrete case it defines the lattice on which 
the canonical Hamiltonian equations are discretized. Being second-order difference equation it 
needs one more initial value (first step of lattice) to state initial-value problem. 

It is interesting to note that the equations (16.61) can be obtained from discrete variational 
equations in Lagrangian framework [TBI El LTSl [JJJ] with the help of discrete Legendre transfor- 
mation 1 201. 



Remark 6.1 Equivalent formulation can be considered for the finite-difference functional 

n 

and a discrete Hamiltonian function ft(£, t + , q + , p). 
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6.2 Invariance of the Hamiltonian action 

Let us consider the functional ( 16.21) on some lattice, given by equation 

n(t,h+,h-,q,p) = 0. (6.7) 



Definition 6.2 We call a discrete Hamiltonian function Ti considered on the mesh (ITl) in- 
variant with respect to a symmetry group generated by the operator ^6. if the action $6.2\) 
considered on the mesh fl£. 7|) is an invariant manifold of a group. 



Theorem 6.3 A Hamiltonian function considered together with the mesh ^6. 7| ) is invariant 
with respect to a group generated by the operator k6. 1\) if and only if the following conditions 
hold 



C D{q l ) +pf D{rf) -X(H)- HD{$ 

+h +h +h 



= o, 

(6.8) 



Xfl(t, h + ,h-, q,p)\ n=0 = 0. 
Proof. The invariance condition follows directly from the action of X on the functional: 

X \Y,PtW + - q l ) - Hh + ) = J2 UtD{q l ) + P tD(v l ) - X(H) - H D(o) h + = 0. 

It should be provided with the invariance of a mesh, which is obtained by the action of symmetry 
operator on the mesh equation (16. 7p . □ 



6.3 Discrete Hamiltonian identity and discrete Noether— type theo- 
rem 

As in the continuous case, the invariance of a discrete Hamiltonian on a specified mesh yields 
first integrals of discrete Hamiltonian equations. 

Lemma 6.4 The following identity is true for any smooth function 7i = T~C(t, t + , q, p + ): 

C+ DW) + Pt nw - x ( n) - n S ({) . « ( S («-) - f - ") 

(6.9) 

^-f(«"+*-f"). 

Proof. The identity can be established by direct calculation. □ 

We call this identity the discrete Hamiltonian identity. It allows us to state the following 
result. 
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Theorem 6.5 The invariant with respect to symmetry operator l}6. 1\) discrete Hamiltonian 
equations 116.6}) possess a first integral 

T = V%-t(n- + h-™ ) (6.10) 

if and only if the Hamiltonian function is invariant with respect to the same symmetry on the 
solutions of equations 116. 6\) . 

Proof. This result is a consequence of the identity (16. 9p . The invariance of the discrete Hamil- 
tonian equations is needed to guarantee the invariance of the mesh, which is defined by these 
equations. □ 

Remark 6.6 Theorem \6.5\ can be generalized for the case of the divergence invariance of the 
Hamiltonian action, i.e. 

C D(q l )+pf D(v l ) - X(H) - HD(0 = D(V), (6.11) 

+h +h +h +h 

where V = V(t, q, p). If this condition holds on the solutions of the discrete Hamiltonian 
equations l}6.6\) . then there is a first integral 

1 = rf Pi - £ (V + ) - V. (6.12) 

Remark 6.7 For discrete Hamiltonian equations with Hamiltonian functions invariant with 
respect to time translations, i.e. 7i = 7i(h + , q, p + ), where h + =t + — t, there is a conservation 
of energy 

oh_ on + 

In this case the discrete Hamiltonian equations 116.6}) are related to symplectic-momentum- 
energy preserving variational integrations introduced for discrete Lagrangian framework in I21\/ . 
Note that 7i is not the discrete energy, it has a meaning of a generating function for discrete 
Hamiltonian flow. 

6.4 Applications 

6.4.1 Discrete harmonic oscillator 

The harmonic oscillator model is very important in physics. A mass at equilibrium under the 
influence of any conservative force behaves as a simple harmonic oscillator (in the limit of small 
motions). Harmonic oscillators are exploited in many manmade devices, such as clocks and 
radio circuits. 

Let us consider the one-dimensional harmonic oscillator 

q = p, p = —q. (6.13) 
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This system of Hamiltonian equations is generated by the Hamiltonian function 

H(t,q,p) = ^(q 2 +p 2 ). 

As a discretization of the equations (16. 13j) we consider the application of the midpoint rule 

q+-q P + P+ P+-P q + q+ 



h, 2 ' h. 



(6.14) 



on a uniform mesh h + = h_ = h. The presented discretization can be rewritten as the following 
system of equations 



h+ \ , , 4 / h+ 



(6.15) 



h + = h_. 

It can be shown that this system is generated by the discrete Hamiltonian function 

2 

K{t,t + ,q,p+) = A _ h2 (q 2 + p\ + h + qp + ). 

Indeed, the first and second equations of (16.61) are exactly the same as those of (I6.15p . The last 
equation of (16.61) takes the form 

2(4 + /^) 16fr + , 2(4 + /£) 2|2 ,, 16fe_ 

-J^W {<1 +P+) -(4^W^ ++ (4^W (g - +P)+ (^W^ = 
Using the first and second equations, we can rewrite it as 

2 +tAtW + p 2 )=0- 



4 + h 2 + 4 + hi 

Therefore, for the case q 2 + p 2 ^ this equation can be taken in an equivalent form 

h + = h- = h. 

The system of difference equations (16.151) admits, in particular, the following symmetries 

d d d d 

Xi = sm(ut)- — h cos(ujt) — , X 2 = cos(ut)- sin(co>t) — , 

Oq op Oq op 

d d" o 1 o 1 o 1 

X 3 = ^r, X 4 = q—+p—, X 5 =p—-q—, 
ot Oq op oq op 

where 

arctan(/i/2) 
" = h/2 ' 

For symmetry operators X\ and X 2 we have the divergence invariance conditions 
C+ D (q) +p + D(v)-X(H)-HD (0 = D (V) 

+h +h +h +h 
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fulfilled on the solutions of equations (16.151) with functions V\ = q cos(ut) and Vi = —q sin (u;i) 
respectively. Therefore, we obtain two corresponding first integrals 

X\ = psm(u>t) — qcos(out), X 2 = pcos(ut) + qsin(ut). (6.16) 

Symmetry operator X3 satisfies the invariance condition 

C+ D(q) +P+ D{r)) - X{H) - H D(0 = 0. 

+h +h +h 



Thus, we get the first integral 



r 4 ( A + h 2 _q 2 _+p 2 4h_ v 



Using the first and second equations of (I6.15p . we can simplify it as 

4 q 2 + p 2 
3 " 4 + hi 2 

Since from the first integrals X\ and X 2 we have conservation 

+ X\ = q 2 + p 2 = const, 
it follows that we can take the third first integrals equivalently as 

X 3 = h_. 

The three first integrals Xi, X 2 , X3 are sufficient for integration of the system (16.141) . We 
obtain the solution 

q = X 2 sin(cL>t) — X\ cos(u;t), p = X\ sin(u;t) + X 2 cos(u;t) 

on the lattice 

U = to + ih, i = 0, ±1, ±2, h = X 3 . 
6.4.2 Modified discrete harmonic oscillator (exact scheme) 

The discrete harmonic oscillator of the preceding example follows the same trajectory as the 
continuous harmonic oscillator, but with a different velocity. This numerical error can be 
corrected by time reparametrization. In this case we will get the exact discretization of the 
harmonic oscillator, i.e. a discretization which gives the exact solution of the underlying ODEs. 
In this case the harmonic oscillator (16.131) is discretized as 

q+-q = q P + p+ p+-p = o i±g+ 
h+ 2 ' h+ 2 (61g) 

h + = h- = h, 

where 

tan(/i/2) 
~ h/2 ' 
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represents a time reparametrization. Similarly to the preceding example it can be shown that 
this discrete model of the harmonic oscillator is generated by the discrete Hamiltonian 

2Q 

H{t,t+,q,p + ) = - - ^ 2h<2 (q 2 +p 2 + + Q 2 h + qp + ). 

The system of difference equations (16.181) admits the following symmetries 

^ . d d d . d 

Ai = smr— — hcost— , X 2 = cost— sint— , 

oq op oq op 

„ 9 d d d d 

X 3 = 7^7, X i = q—+p— 1 X 5 =p—-q—. 
ot oq op oq op 

For symmetries X\ and X 2 , which satisfy the divergence invariance condition (16.111) with func- 
tions V\ = qcost and V 2 — —qsint, we obtain two first integrals 

Xi = psint — q cost, X 2 = pcost + qsint. (6.19) 

The operator X 3 satisfies the invariance condition (16.81) and provides us with the first integral 
X3, which (similarly to the preceding example) can be taken in an equivalent form 

J 3 = h- (6.20) 

The scheme ( I6.18P gives the exact solution of the harmonic oscillator, which can be found 
with the help of first integrals X\ and X 2 as 

q = I2 sin t — X\ cos t, p = X\ sin t + X 2 cos t. 

This discrete solution is given on the lattice 

U—to + ih, i — 0, ±1, ±2, h = X 3 . 

The exact schemes for two- and four- dimensional harmonic oscillators were used in [22] to 
construct exact schemes for two- and three-dimensional Kepler motion respectively. 

6.4.3 A non-linear motion 

We consider a difference analog of equations (15.111) . and choose 

H(t,t + ,q,p + ) = ~(pl + (6.21) 
Then, in accordance with (16. 6p we obtain the discrete Hamiltonian equations: 

D(q)=p+, D(p) = ^, 
+h +h q 6 

(6.22) 

2 1 2 1 

P+ + -o = P + ~ 
q 2 qi 
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It is easy to check invariance conditions for Tt with respect to symmetry operators X\ and 
X 2 , given in (15.121) . Application of Theorem 16.51 for these symmetries yields first integrals 




) 



I 2 = qp-tip 2 + — 



) 



(6.23) 



Therefore, the solution of the discrete system (16.221) satisfies the relation 



I2 = qp + Ztlx 



in all points of the lattice. 



7 Conclusion 



The goal of the present paper is to present a method to find first integrals of canonical Hamil- 
tonian equations and to establish a way to preserve Hamiltonian structure in finite-difference 
schemes. To achieve this we use invariance of the Hamiltonian action functional and its rela- 
tion to first integrals of canonical Hamiltonian equations. The conservation properties of the 
canonical Hamiltonian equations are based on the newly written identity (called the Hamil- 
tonian identity). This identity can be viewed as a "translation" of the well-known Noether 
identity into the Hamiltonian framework. The identity makes it possible to establish one-to- 
one correspondence between invariance of the Hamiltonian and first integrals of the canonical 
Hamiltonian equations (the strong version of Noether's theorem). 

The variational consequences of the Hamiltonian identity make it possible to establish nec- 
essary and sufficient conditions for the canonical Hamiltonian equations to be invariant. These 
conditions make it clear why not each symmetry of the Hamiltonian equations provides a first 
integral. 

The Hamiltonian version of Noether's theorem, formulated in the paper, gives a constructive 
way to find first integrals of the canonical Hamiltonian equations once their symmetries are 
known. This simple method does not require integration as it was illustrated by a number of 
examples. In particular, we considered equations of Kepler motion in various dimensions. The 
presented approach gives a possibility to consider canonical Hamiltonian equations and find 
their first integrals without exploiting the relationship to the Lagrangian formulation (see, for 
example, [23]). 

The approach developed for the continuous case was applied to discrete Hamiltonian equa- 
tions, which can be obtained by variational principle from finite-difference functionals. Sim- 
ilarly to the continuous case we related invariance of discrete Hamiltonian functions to first 
integrals of the discrete Hamiltonian equations. In particular, energy conserving numerical 
schemes can be obtained as discrete Hamiltonian equations generated by Hamiltonian func- 
tions invariant with respect to time translations. 

The results presented in the paper can be used to find first integrals of continuous and 
discrete canonical Hamiltonian equations. They also provide guidelines how to construct con- 
servative finite-difference schemes in Hamiltonian framework that is important in numerical 
implementation. 
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